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The current supernova detection technique used in IceCube relies on the sudden deviation of the 
summed photomultiplier noise rate from its nominal value during the neutrino burst, making IceCube 
a ~ 3 Megaton effective detection volume - class supernova detector. While galactic supernovae can 
be resolved with this technique, the supernova neutrino emission spectrum remains unconstrained 
and thus presents a limited potential for the topics related to supernova core collapse models. 

The paper elaborates analytically on the capabilities of IceCube to detect supernovae through the 
analysis of hits in the detector correlated in space and time. These arise from supernova neutrinos 
interacting in the instrumented detector volume along single strings. Although the effective detection 
volume for such coincidental hits is much smaller (> 35kton, about the scale of SuperK), a wealth 
of information is obtained due to the comparatively low coincidental noise rate. We demonstrate 
that a neutrino flux from a core collapse supernova will produce a signature enabling the resolution 
of rough spectral features and, in the case of a strong signal, providing indication on its location. 

We further discuss the enhanced potential of a rather modest detector extension, a denser array in 
the center of IceCube, within our one dimensional analytic calculation framework. Such an extension 
would enable the exploration of the neutrino sky above a few GeV and the detection of supernovae 
up to a few lOO's of kilo parsec. However, a 3 — 4Mpc detection distance, necessary for routine 
supernova detection, demands a significant increase of the effective detection volume and can be 
obtained only with a more ambitious instrument, particularly the boosting of sensor parameters 
such as the quantum efficiency and light collection area. 



I. INTRODUCTION 

Much of the underlying physics of a core collapse su- 
pernova remains unexplored, the mechanism leading to 
the explosion not well understood. Though it is clear 
that neutrinos play a major role in the dynamics of the 
explosion, carrying away 99% of the energy emitted by 
the dying star. The observation of neutrinos in the prove- 
nance of SN 1987A in several experiments [lH3] has con- 
firmed the crucial role they play in the overall picture. 
The limited number of detected neutrinos, however, has 
not provided enough statistics to constrain current super- 
nova models. This may only be achieved with a detector 
performing routine SN detection, about one per year. 

The idea of using IceCube as a supernova detector was 
first discussed in [H [S]. The method is based on the 
sudden increase in the photomultiplier count rate on a 
timescale of the order of 10 seconds: qualitatively, the re- 
lease by a supernova at d = 10 kpc with a total amount of 
energy of 10^^ MeV in neutrinos of about 20 MeV corre- 
sponds to a fiux of 7 • 10^** m~^ In its completed configu- 
ration and with its data acquisition system, IceCube has 
about 3Mton effective detection volume and a reach of 
about 50 kpc at ha C.L. detection level. Here, we extend 
that method and assess the potential of IceCube and sev- 
eral prospective detector configurations by considering 
multiple hit detection in a single or neighboring sensors 
from individual neutrino interactions. We demonstrate 
that multi-hit modes enable the extraction of new ob- 
servables related to the features of the neutrino emission 
spectrum and the supernova location, not accessible by 
means of the original method. This new data stream 
could help study questions related to supernova dynam- 
ics (neutronization burst, accretion and cooling phase. 



see [7] and references therein) or fundamental questions 
of particle physics, related to the nature of the neutrino 
[H in]. The observation of electron neutrinos from the 
early neutronization burst could disentangle the neutrino 
hierarchy provided a sufficiently large mixing angle 
6*13 as indicated recently [TUj . 

These new detection modes are conceptually situated 
between a full reconstruction of the events and the cur- 
rent search method balancing the simplicity of the latter 
with the richness in information of the former. We show 
that the loss of information in the current method is not 
necessary and that a very elementary reconstruction, im- 
plemented through the coincident hit methodology, can 
enhance the information obtained during a supernova ex- 
plosion. In fact, it would appear natural to extend the 
new detection modes to a full reconstruction of super- 
nova events in the case of a very densely instrumented 
volume. 

According to [TH [13] , the estimated SN rate from our 
galaxy is about 0.03/yr. It increases to 0.1 /yr within 
1 Mpc and about 1/yr within 4Mpc. We will see that 
extending the reach of a future detector to such a large 
distance is quite challenging given the quick decrease of 
the interaction density with distance. SuperK, with its 
32 kton fiducial volume, remains at the moment an in- 
strument of choice jU] . IceCube jT^ with its 86 strings 
and this new detection methodology could become a se- 
rious competitor, even surpass it, provided a deep and 
denser core extension installation. There are several good 
reasons for such an extension beyond an enhanced poten- 
tial for SN detection, mainly the detection of individual 
high energy neutrinos down to several GeV related to 
multiple physics aspects, such as neutrino sources in the 
Southern sky [TB] [T7] , indirect detection of dark matter 
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down to lower mass [TS^EU] and neutrino particle physics 

This contrasts with the potential of other projected 
lOOkton-class detectors, e.g. liquid scintillator LENA [55] 
or liquid Argon detectors [23j . which will have an ex- 
tended potential for SN detection but comparatively 
small fiducial volumes at energies beyond GeV. Espe- 
cially the latter are complementary to water or ice-based 
Cherenkov detectors, since they are more sensitive to the 
electron neutrino channel, providing a better access to 
the physics during the neutronization burst [24] . Com- 
bining the measurements from these detectors, which dif- 
fer not only in their sensitive material, but in their geo- 
graphic location, would enable the study of matter effects 
and neutrino hierarchy by observing neutrinos from the 
same supernova from different viewpoints with a high de- 
tection probability. 

Finally, a large detector will not only increase the dis- 
tance reach for SN but also provide large event statistics 
for galactic SN with possible resolution for their loca- 
tion, leading through a SN alert system (SNEWS f5S^) 
to the optical detection of obscured and otherwise missed 
galactic SNe. 

In Section 2, we define the studied configuration and 
the benchmark supernova neutrino flux. We present the 
effective volume analytical calculations in Section 3 and 
their consequences in Section 4. 



II. DEFINITIONS AND METHODOLOGY 



Detectors and sensors 



Configurations 

In the following, we consider several detector configura- 
tions [55] : 

a. The baseline IceCube configuration (IC), with a 
number of sensors iVoM — 4800, arranged i = 17 m apart 
along 80 strings. 

b. The Deep Core configuration (DC), with iVoM — 
320, L — Ixn apart on 8 strings. 

c. The Dense and Deep Core configuration: An 
imaginative futuristic low energy extension in alterna- 
tive configurations denoted DDC, DDC47r, DDCY^ com- 
posed of 24 strings with 150 sensors each, iVoM = 3600, 
and with a spacing between modules of i = 3 m. In 
configurations with Ati subscript, the sensitive area of 
the spherical module has a 47r acceptance, for instance 
the modules designed for KM3NET [27]. The DDC^.^ 
configuration has also a larger sensitive area and higher 
Q.E. In addition, we assess the potential of the full accep- 
tance configurations with extremely low noise modules 
(100 Hz). 

A short local coincidence is the occurrence of two or 
more hits within specific time windows 5Vs in a single 
sensor or in two neighboring sensors, also referred to as 



optical modules. The systematic recording of these co- 
incidences would require the implementation of a dedi- 
cated data acquisition mode and its corresponding data 
stream. We study in this paper several detection modes 
* ^ 2 and j < 1, referring to the coincidental detec- 
tion of i hits with a module and j hits with the module 
just above along a sensor string. 

In this paper we do not consider coincident hits be- 
tween strings, so string spacing is not a parameter in 
this study. We refer therefore to our study as ID ana- 
lytic approximation. We also restrict the study to nearest 
neighbor (NN) coincidences. 

Sensor geometry 

The geometry of the optical modules for the three de- 
tector configurations IC, DC and DDC is illustrated in 
Fig. [T] The modules [28] are composed of a 0.5 inch en- 
veloping glass sphere (to withstand the pressures at the 
large deployment depths, in particular during the ice re- 
freezing process) with an external radius R = i?oM = 
13 inch, enclosing a photomultiplier R7081-02 made by 
Hamamatsu Photonics 33 . The glass and contact gel 
between the PMT and the glass introduce a 350 nm lower 
cutoff. 

In all configurations except DDC47r, the photocath- 
ode sensitive surface of the sensor covers a fraction of 
the sphere slightly smaller than the half sphere. The re- 
sulting cutoff angle on the sensitive surface is given by 
^max — '?i'/2, where 0,nax is the limit on the incoming di- 
rection of Cherenkov photons at which they can intersect 
the sensitive surface of the module. Given the specifica- 
tions of the IceCube photomultiplier, we can estimate 

TT ^max ~ 25 . 

Viewed under an angle 0, the projected area of the 
sensitive surface is 



A = AoM{e) = nRR{e), 



where 



R[e) = R- 



1 — cos ( 



Hi 



(1) 



0) (2) 



for a module sphere looking down and with a photocath- 
ode sensitive area cutoff angle ^max ^ 7r/2, as illustrated 
on Fig. [l] 

Sensor photon detection efficiency 

The photon detection efficiency (PDE) of the sensitive 
surface x(-^) has been measured somewhere |50| . it com- 
bines the PMT Q.E., the transparency of the surrounding 
glass and the contact gel which couples the glass to the 
PMT. Averaged over a Cherenkov spectrum between 
300 and 650 nm, and accounting for an amplitude trigger 
threshold keeping ethr =90% of the p.e.'s, it yields an 
integrated PDE 



ethr 



x(A)/A2 d\ ^ 7.3% 



(3) 
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FIG. 1: Illustration of the sensor geometry and definition of 
6 and 0^3.^ . 



for the regular IC sensors. The high quantum efRciency 
sensors deployed in DC are assumed to have a PDE 
30% higher [26]. The same assumption is made for the 
DDC and DDC4^ configurations. In contrast to previ- 
ous works [5j, (x) is significantly lower, due to a regular 
optical module peak PDE of about 14% at 420 nm. 

Although {x) is a fair approximation for our analyti- 
cal estimates presented in the next section, a toy Monte 
Carlo simulation should implement the complete wave- 
length dependence of x(-^) given that the propagation 
of photons in the ice medium is described in terms of 
wavelength-dependent scattering and absorption lengths 
-'^scatt and Aabs 133 ■ This dependence in the propagation 
cannot be disentangled from that of the PDE. 

Sources of background noise 

The following types of background may limit the appli- 
cability of the coincident hit detection methodology: 

a. Cosmic-ray induced muons. A fraction will travel 
down to the detector and produce coincident hits. We as- 
sume that most of these events are muons dying in the 
upper layers of the detector: they barely enter the fidu- 
cial volume and produce light on a very limited number of 
external sensors. We assume that this event class can be 
adequately rejected. Some higher energy muons however 
penetrate deeper into the detector and although most of 
them can be vetoed given their specific topology, some 
may leave a signature well inside the instrumented vol- 
ume mimicking the coincidence signature under investi- 
gation. We would expect the resulting rate to be minimal 
and thus neglect it, but this specific topology should still 
be further studied separately by means of the experimen- 
tal data, to properly substantiate this statement. 

b. Low E atmospheric neutrino-induced muons or 
electronic showers. This background is constituted by 
low energy neutrinos in the 10 MeV — 1 GeV range [551 - 



157] . The interaction rate is very limited, well below 1 Hz, 
in a Gton ice volume. 

c. Sensor correlated noises at fis and ms time scales, 
due respectively to photomultiplier afterpulse and ra- 
dioactive decay chains |45j. The afterpulse components 
of the dark noise come with well defined [JH] |46] delays 
after the original photo-induced signal, larger than the 
time windows used below for coincidental hit detection 
within a single sensor, so that (1) it does not pollute and 
therefore impede the coincidental hit detection method 
for this specific channel, (2) it can be filtered out by 
means of an offline analysis, leading to an effective noise 
rate significantly smaller than measured (dark count re- 
duced from 580 Hz to r = = 440 Hz 28, 29 ). This 
reduced dark count rate value is thus considered to calcu- 
late the probability of the noise. Note that the deadtime 
introduced by this afterpulse filter is negligible (< 0.5%). 

The effect of ms timescales of correlated noise is more 
complex and it is not clear whether it would affect the po- 
tential of the new methodology when applied to specific 
channels with multiple hits in a single sensor. 

We define the hit times {ti}i, a hit time window by its 
width {St)ij and an arbitrary time shift t^ . Hits i and j 
belong to the time window if < {5t)ij/2. We 

restrict the study to several qualitatively complementary 
NN detection modes, 

• single sensor hit: mode 10 or "1-1-0"; 

• two photons detected in the same sensor: mode 20 
or "2-HO"; 

• a single hit detected in two NN sensor: mode 11 or 
"1+1"; 

• three photons detected at times ti, t2 (in the same 
sensor, ti < t2) and ^3 (in the sensor above): mode 
21 or "2+1". 

The noise rates of two hit coincidences and r^^ 



given by 



rii=r2(,Si)i3andr20^2r2(,5t)i2 (4) 



(the latter factor two because hits at a single module 
are indistinguishable). The coincidental noise rate with 
two hits in the lower sensor and a single hit in its upper 
neighbor is 



r''^2r' ((5<)i2 {5t) 



13- 



(5) 



Note that in the IceCube detector, only the lower hemi- 
sphere of the optical module is sensitive and therefore, a 
hit in the upper sensor coming after a hit in the lower 
sensor is strongly preferred (i.e. an optimal time win- 
dow will have iis > 0). This is not true for some of the 
futuristic detector configurations which have a uniform 
sensitive surface: in this case, we will choose ^13 = 0. 

The detector coincidental noise rate scales linearly with 
the number of sensors (the number of independent pairs 
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is almost equivalent to the number of sensors, provided 
that the number of sensor along a string is ^ 1). 

IceCube has a system that records local coincidences 
within a 2/is (±1 /xs) time window up to the next-to- 
nearest neighbor (there are four such neighbors) [SU] , 
The rate of such coincidences, called hard local coin- 
cidences (HLC), should be therefore of the order of 
r^^^ — 4 X 2 X ri^s w 1.55 Hz per optical module. 
For our purpose, we consider the rate of coincidental 
events within a much shorter coincidence time window 
with At = 30 — 150 ns involving one neighboring sensor 
corresponding to a coincidental noise rate of r-so- 150 ns ~ 
0.006 — 0.03 Hz per sensor. 

A SN signal becomes significant when the number of 
coincidental events integrated over the typical duration of 
the SN neutrino emission (about 3 seconds) takes values 
larger than the fluctuations of the noise rate r'-' , corre- 
sponding to the various event classes. 

B. Supernova neutrino emission models 



is interesting to note, as illustrated in Fig.[2j that the pa- 
rameters a and {E) can be very well fitted with Eq. [t] in 
order to match the integrated neutrino flux and positron 
density up to a certain time. 
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We have chosen to study several flux models of rather 
modest luminosity. Our main benchmark flux model is 
the "Sf model (also referred to as "Garching" model), 
an 8.8 solar mass electron capture supernova, discussed 
in [S| and which may represent a large fraction of all 
supernovae. 

Furthermore, in this paper we will only consider the 
predicted anti-electron neutrino flux, focusing solely on 
the inverse beta decay (IBD) reaction in ice, i.e. I'c+P — > 
n-t-e"*". For one, in water or ice, which are natural media 
for high energy neutrino telescopes, the (anti-)electron 
neutrino interaction with oxygen has a much higher en- 
ergy threshold and a quite significantly smaller CC cross 
section at the energies of interest [31] . This largely makes 
this channel subdominant for the interaction rate in the 
detector compared to the IBD reaction yield. Moreover, 
a large fraction of the neutrino energy is absorbed into 
binding energy, reducing the Cherenkov light yield and 
dramatically suppressing the potential of this reaction 
channel given that our methodology relies on the co- 
incident detection of Cherenkov photons in neighboring 
sensors. We also neglect the elastic interaction channel 
j/c + e^ — !■ Pc + which has an even lower cross section 
than the oxygen channel |32) . 

The anti-electron neutrino flux luminosity is given by 



d$ 



where 



1 



m 

And^ Jo {E,){t) 



fa(t),{E^)(t)dt (6) 



fc 



(l + a)i+" / E^ y (1+^)^ 

^<-">(^'')=(£;.)r(i + a) [W)) ' 



and with the time-dependent model parameters luminos- 
ity L{t), average neutrino energy {E^){t) and pinch pa- 
rameter a(t). Their time evolution can be found in [6 . It 



FIG. 2: Positron spectrum, integrated up to t. Dashed line is 
the result of fitting a = 4.77 and (Ec) = 16.97 MeV. 

The neutrino flux integration time t, which closely op- 
timizes the potential in terms of distance reach of the co- 
incident hit method is derived by maximizing the signal 
to the square root of noise ratio, assuming an exponen- 
tial luminosity decrease, in good qualitative agreement 
yt > Is, with tQ « 2.35s, of the luminosity evolution 
in the Garching model (a usual luminosity decrease in 
{1 + t/toy^ does not fit wefl). 

The argument could be more sophisticated, e.g. by us- 
ing a refined likelihood temporal analysis. 

The integrated signal (the supernova neutrinos pro- 
ducing an event) up to time i is S{i) — soe^*/*°dt — 

soto{l - e-*/*»). The detector noise is assumed gaussian 
with count fiuctuation of a complete detector configura- 
tion y/Jl = \J NoMfi within an integration time i, where 
r is the noise rate of a specific mode, i is determined 



leading to 



t ■ 



d S{i) 
di JTl 



|(e*7*o 







!)• 



(8) 



(9) 



We extract t — 2.95 s for Iq — 2.35 s and we note that t 
is independent of sq (proportional to the luminosity), of 
the noise rate r and of the significance of the observation. 

The luminosity at early times is larger than in this 
description, so that this integration time t constitutes 
a conservative upper bound regarding the potential of 
the method: if more precisely assessed, t and thus the 
integrated noise would be reduced. 
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The integrated anti-electron neutrino luminosity 
within this time-span t is L = 2.13 • 10^^ erg, and the 
flux = 0.89 • 10^^ m^^. The average neutrino energy 
is (£^p^) = 12.5 MeV. L is about three quarter of the to- 
tal luminosity. It is important to note, when comparing 
the limiting distance up to which a SN can be discov- 
ered, that the total luminosity is about 56% of the one 
considered in [TB] . 

We have also calculated the outcome of the novel de- 
tection methodology for neutrino spectra which obey a 
Fermi-Dirac distribution, 

/,,T(i?.) = iV(,7,T)^^^^^ (10) 

where N{t],T) = -{2TLi3{-e'')y\ We have normal- 
ized the amount of energy released over t to be equal 
to the Sf model, for two temperature, T — 5 MeV 
and T = 6.5 MeV. Although the average neutrino en- 
ergies, (£^pJt=5MgV = 15.8 MeV and (i^pjT=6.5MoV = 
20.5 MeV, are larger resulting in smaller flux than in the 
Sf model, the interaction density is increased by more 
than 50% due the IBD cross section increase turning into 
a considerable increase in the coincidental hit rate (by a 
factor > 2) as will be seen. 

C. Discovery potential and sensitivity 

We assess the discovery and learning potential by re- 
questing the signal over square root of noise ratio s = 5 
and s = 1 for the C.L. in a unit, respectively, using the 
Poissonian distribution p^: given a mean number of ex- 
pected coincidental noise events //, a minimum number 
of > TO events is required to reach significance > s, 

m 

min{m| p^(0) + ^p^(fc)>erf(^)}. (11) 
" k=l 
The distance reach is therefore given by 

c'lokpc = y/ns/im- fi), (12) 

where ^lokpc distance in 10 kpc units and is the 

number of signal events from a SN explosion at 10 kpc. 

An advantage of the detection method relying on co- 
incidental hits arises from a more gaussian distributed 
coincidental noise rate (while the single hit noise rate dis- 
tribution has non gaussian tails, requesting coincidences 
dilutes them). The false alarm rate (assuming gaussian- 
ity) for s = 5 is roughly f/erfc(5/\/2)/2 — 3yr^^. 

s = 1 is a requirement too weak and not adequate for 
IceCube to trigger as a standalone detector, but where 
IceCube may help establishing stronger constraints by 
enhancing the quality of an observation: it sets the dis- 
tance scale for IceCube to contribute in supplementing 
information to observations made by other instruments 
detecting e.g. a neutrino or an optical signal. In the 
latter case, of course, IceCube would have to buffer the 
relevant information because of the few hours delay be- 
tween the neutrino and the optical signal. 



III. ESTIMATE OF SINGLE AND COINCIDENT 
HIT DETECTION RATE 

The positron effective detection volume as a function of 
the neutrino energy is the product of the number of target 
protons in the interaction volume and of the fraction of 
uniformly distributed positrons leading to a detection. It 
is an important and well defined concept for our analysis 
because the positron tracks are nearly point-like with a 
path length of about 10 cm in average, a small distance 
compared to the mean distance between a positron and 
a sensor in which it leaves a signature. This enables the 
estimation of the response of the detector to a given flux 
of low energy neutrinos. 

The positron effective volume for coincidence of order 
fc at a single module is denoted V^{E), where k = 1,2 
for the detection of a single hit resp. of two coincident 
hits. We write V^f^ (E) for a detection of order k in one 
sensor and of order k' at the NN. 

The effective volume depends only slightly on the in- 
coming neutrino direction due to the proximity of the 
average neutrino energy with the energy at which the 
anisotropy in the positron emission (transition between 
backward and forward scattering) from the IBD interac- 
tion vanishes. Note however that selecting events in the 
high energy tail of the neutrino emission, where forward 
scattering is favored, may grant access to the direction 
of a nearby supernova. This can be done to some extent 
by means of the coincident hit time analysis. 

The current supernova detection strategy with Ice- 
Cube relies only on V^ff{E). We analyze here the power 
of the new strategy by relying more generally on V^g (E) . 
We emphasize as well the importance of certain ratios, 
such as V;|(£;)/V;ij(£;) or V^iiE)/V^g{E), which open 
new perspectives and grant access to new observables, 
e.g. the average positron energy. 

For a given supernova neutrino emission model M, 
time-integrated up to time i, the effective volume defi- 
nition can be extended to embed the spectral shape of 
the positron energy distribution dpj^ /dE, 

where 

P^{d)= J -jf{E,d)dE (14) 

and which follows from the definition 

N^^>''^'{d) = p^{d)Vj^''"''. (15) 

From now on, we drop the Ai or replace it by one of the 
specific models. The presence of the tilde denoting time- 
integration necessarily refers to a model and confusion in 
the notation is avoided. 

The differential positron density dp^+ / dE^ is extracted 
from the incoming neutrino flux d^^/dE,^ by folding it 



6 



with the differential IBD cross section d(j/dE^+ given by 
Eq. (10) in ^ (accurate up to NNLO in e), 

(16) 

and the total number density of positron interactions 
(same luminosity for all three models), 

0.074 m-3, Sf, 
Pe+(10kpc) = { 0.096 m-3^ Tfd = 5MeV, (17) 
0.12 m-3, TFD = 6.5MeV 

The proton target density is given by pp = 2 piccNa/Micc, 
where p\cc — 0.91 g/cm^, Mice = 18.015 g/mol, and Na is 
the Avogadro constant. 

The number of expected ev ents belonging to the var- 
ious classes 7V''+ (see Eq. 15), obtained by integration 
of 
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FIG. 3; Number of emitted photon / 10 / MeV by a positron 
and its byproducts between 300 and 650 nm (Geant4) (solid 
line) and positron range parametrization above Cherenkov 
emission threshold in mm/MeV (dashed line) w.r.t. the 
positron energy. 



dA^, 



fcfc' 



-{E,, d), 



(18) 



B. Effective detection volume 



also defines the relevant observables rj.-'f^, = N^^ /N-'l , 
some of them being particularly sensitive to the average 
energy of the neutrino spectrum. 



A. Cherenkov emission 

In order to proceed with the effective volume estimates, 
we need to approximate the average number of detectable 
Cherenkov photons emitted along the track of a positron 
with energy E, which could naively be done by combining 
the Franck-Tamm [32] (/? = 1), the average photomulti- 
plier Q.E. and positron range formulas. 



N(E) = 2TTaR{E) 1 



1 



A2 



(19) 



where the integrals for average Q.E and Franck-Tamm 
formula run between 300 and 650 nm. R{E) is the 
positron track length above Cherenkov emission thresh- 
old and given by 



(20) 



R{E) = In 2 In (1 + E/{E^ In 2)) A^'^^^Vpi, 



where 



36.08 g/cm2, piee = 0.92g/cm2 and E^ 



61GMeV/(Z-|- 1.24) = 72.1 MeV gS]. 

E.g. at £; = 10 MeV, R{E) w 5.5 cm and the aver- 
age number of detectable photons (IC) by a positron is 
A^ K, 135. However, Eq. [19] does not take into account 
Cherenkov photons which are emitted by secondary par- 
ticles, e.g. knock-off delta electron and positron annihi- 
lation. Results from a Geant4 simulation show that the 
total number of emitted photons remains almost constant 
(c.f. Fig. |3]). In the following we will use the Geant4 pho- 
ton yield results for the calculations. 



The fraction e(r, 9) of the 47r sr solid angle occupied by 
the sensitive surface of the module, viewed under angle 
6 and at distance r is 



1 



1 



e r, 



Vi + RR{e)/r^ 



(21) 



(we have neglected the spherical shape of the sensor but 
this remains an excellent approximation). 
When r ^ i?, this expression reduces to 



A 



DOM 

47rr2 



(22) 



The expectation number of detected photons is 

p{E,r,e) = N{E)e{T,e) (23) 

We have now all the necessary ingredients to express 
the coincident hit effective volumes, 

V^jf{E,) = 2ti / r^Ar / dcos(?/(fc, r, 9, E,). (24) 

fik, r, 0, E,) = J (^1 - J2PKr,e,E^)il)^ (25) 

is the density of probability of recording > k hits from a 
positron with energy E^, where e, e^) {k) is a Poisson 
distribution with mean value 

/i(r, 0, E,) = AN{EMr, &) exp (-r/Aabs) (26) 

For the effective volume for a sensor with k hits and 
NN sensor with k' hits, we can write the more generic 
expression. 
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V^ff{E) = 2tt r^dr / dcos6' 

X if{k,r,0,E,)f{k',r',e',E,) (27) 

where r' = y/r'^ + L"^ + 2Lr cos 9 and 008 6*' = [L + 
rcos9)/r'. 

According to [5 and for the purpose of this esti- 
mate, we have made the quahtative assumption that the 
Cherenkov photons are emitted in correlated directions, 
roughly over Trsr, or a quarter of the sphere, the approxi- 
mate size of the Cherenkov cone. We therefore effectively 
implement the anisotropic emission with the inclusion of 
several factors 4. In order to obtain simple expressions, 
we have neglected the spacing between the modules which 
do affect the probability of two NN sensors to be in the 
pool of photons. This may have a significant impact on 
the effective volume. We qualitatively argue that most 
interactions leading to a coincidental hit detection occur 
slightly below the lower sensor, i.e. 9 is in general close 
to 9'. However, this argument does not hold for An sen- 
sitive modules. In this case, the expressions derived here 
are quite approximative. 

The correlated probability of photon detection is jus- 
tified at short distances from the sensor, as the photons 
are not yet diffused. This correlation gradually vanishes 
at larger distances and fully disappears at an effective 
scattering length distance. Fortunately, our assumption 
is likely quite accurate for NN coincidental hit detection 
from low energy events, because the detection probabil- 
ity of > 1 photon is sizable only at short distances and 
the photon directions are not yet decorrelated (i.e. the 
fraction of contributed events beyond a relatively short 
distance is negligible, due to a negligible fraction of the 
subtended areas of the sensors). 

The validity of the effective volume formula for sin- 
gle hits with the introduced anisotropic emission holds 
for arbitrary distance: although large values of the inte- 
gration variable r (i.e at distances where the Cherenkov 
photons and the initial positron directions are decorre- 
lated) contribute significantly to the integral, there is an 
equivalence at large distances between isotropic emission 
- uniform low photon expectation value and anisotropic 
emission - higher photon expectation value within the 
TT sr cone. Stated differently, the equivalence at large dis- 
tances (i.e. /i small) f{fi) « 4/(/i/4) (see the expression 
for /(^), Eq.[25]). 

In Eq. |26[ we have introduced a damping factor related 
to the absorption length. One may question whether 
we should consider instead the attenuation length as ap- 
propriate. This can be understood from first principles 
(phase space conservation) : photons do not propagate to 
such large distances as they scatter and will only reach 
an average distance of about an attenuation length from 
the sensor. However, the subtended sensor area for the 
photon detection decreases with distance and the large 
absorption length considered compensates for the shorter 



effective distance reached by an average photon but see- 
ing a larger subtended photo-sensor area. 

The results for the effective volume in the various de- 
tection modes are shown in Fig. |4j It is noticeable that 
while the ratio of the effective volume to the positron en- 
ergy for single hit detection mode behaves rather flat and 
even eventually decreases at high energy, it continually 
increases for coincidental hit detection modes. The rate 
of the increase is fastest for the 2 -|- 1 mode and faster for 

1 + I than 2-1-0. This behavior reflects in a correspond- 
ing increase of the average positron energy for the chosen 
neutrino emission model. For the four different topolo- 
gies 1 -I- 0, 2 -I- 0, 1 -|- 1, 2 -|- 1, the relative variation rate 
with energy is constant, being independent of the spe- 
cific detector configuration, the absolute scaling of the 
effective volume, which in turn depends on the quantum 
efficiency, and the sensor coverage and spacing. 

The effective detection volumes in the DDC configura- 
tion are significantly higher when the sensor is sensitive 
over 47r, a factor about three for 1 + mode, four for 

2 + mode and six for 1 + 1 and 2 + 1 modes around the 
Garching energy {EjyJ- 

Fig. [5] shows the detection probability w.r.t. the IBD 
neutrino interaction location for DC and DDC47r config- 
urations. 

It is clear from these formulas that, while we constrain 
ourself to a limited number of detection modes, it is triv- 
ial to extend the analysis in order to obtain effective 
volumes for arbitrary coincidental modes kk', including 
those beyond our one dimensional analysis (i.e. modes 
with hits on different strings which become relevant for 
a hypothetical configuration with small string spacing). 



IV. RESULTS 
A. SN detection distance reach 

With the effective volumes estimated above, we now 
have all the ingredients to calculate the number of single 
and coincidental hit detections. In Table |T] we provide 
the number of detections and the distance reach for unit 
signal over background information and for a 5a C.L. de- 
tection, calculated for the various configurations with the 
Sf model, and a comparison of effective volumes in Ta- 
ble [nj Note that in the case of the Fermi-Dirac spectrum 
with T = 6.5 MeV and the luminosity Lbk calculated by 
[l3] , and large area low noise sensors, the 2 + 1 detection 
mode reaches = 2.9 Mpc and d^^ = 7.2 Mpc. 

The description of the distribution of hit arrival time 
differences (see [17| and references therein) accounts for 
the characteristics of the ice, which influences the pho- 
ton propagation. The shape of this distribution depends 
on the distance of the sensor from the photon emission 
location as illustrated in Fig. [6j which shows the arrival 
time difference of hits at one sensor and its (upper) neigh- 
bor (non uniform sensor). We notice that the shape in 
2 + mode depends weakly on the detector configura- 
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FIG. 4: Upper / middle /lower plots: effective volumes for sin- 
gle hit / coincidental hit / next order coincidental hit. Dotted, 
dashed, dot-dashed and solid curves resp. for DDC4,r, DDC, 
DC and IC. Middle plot, thick / thin lines for resp. 2 -f and 
1-1-1 modes (note that DC and DDC results are the same for 
1 + and 2 + modes). 



tion, and that a trigger based on the difference of arrival 
time of hits within bt — 50 ns catches most of the events 
(>95%). We use this value for all detector configura- 
tions. For the 1 + 1 mode, we optimized the cut on the 
time difference window bt — tj, ~ t\ (^1,3 are resp. the 
arrival times at the lower and upper sensor) in terms of 
d^'^, and rescaled ris by the fraction of events surviving 
the time cuts bt = 50, 100, 150 ns, resp. for DDC's, DC 
and IC detector configurations, accordingly. 

As we pointed out earlier, the distance reach with 




-30 \ 1 1 1 1 7 

5 10 15 200 5 10 15 200 5 10 15 20 

Rxv iml 

FIG. 5: Detection probability as a function of the interaction 
for 1+0 and 2+0 modes in DC and 1 + 1 mode in DDC4,r. The 
log of the equiprobability contours correspond to the values 
explicitly shown on the right plot. 

DDCyj" and sensor noise reduced to about 100 Hz has 
a reach larger than 2 Mpc in the 2 + 1 mode, using the 
benchmark luminosity of [13], twice the IC Q.E., and 
modules with 50% larger radius. It is eventually lim- 
ited by the signal statistics and not by the noise, so 
that it constitutes an upper limit of this > 0.5Mton 
(2 + 1 mode) configuration. In order to perform rou- 
tine SN detection, it is necessary to have > 4 Mpc. 
The detector extension must therefore be significantly 
more ambitious and provide at least an effective volume 
of 2Mton, assuming that the noise is not limiting the 
reach. For an IceCube extension that would mean the 
installation of > lO'OOO highly efficient modules, proba- 
bly with a very small string spacing, in order to collect 
a substantial amount of events with hits spreading over 
more than one string. Reducing further the spacing be- 
tween sensors slightly improves the reach, e.g. to 2.2 Mpc 
with 2 m, 3 m being quite close to optimal w.r.t. the stud- 
ied signatures. In conclusion, the realistic implementa- 
tion of a supernova detector with a reach up to 4 Mpc 
seems quite demanding from the point of view of costs, 
deployment time, complexity and R&D, but not unmoti- 
vated given the great physics return. Not to mention the 
improvements it would bring for neutrino studies in the 
1 — 100 GeV range (oscillations, indirect search for dark 
matter, southern sky and low energy neutrino sources). 

In conclusion, other projects such as a 0.5 Mton Hyper- 
Kamiokande [3H]j or the aforementioned lOOkton LENA 
[22] would have a quite limited potential as well w.r.t. 
the detection and characterization of spectral neutrino 
light curves from extra-galactic supernovae. We there- 
fore confirm that a large multi-Mton effective volume is 
necessary to this purpose, such as the 5 Mton DEEP 
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conf. 


1 7 ' 

kk 


Vcff/m 


Us 


ribg 


a /kpc 




1+0 


764 


2.70 lO'' 


6.210" 


46.5, 104 


IC 


2+0 


7.4 


2.085 10^ 


274 


49.9, 108 


1+1 


1.7 


450 


411 


20.9, 45.4 




2+1 


0.12 


26.5 


0.02 


29.8, 52.0 




1+0 


990 


2.3410" 


4.1510" 26.9, 60.1 


DC 


2+0 


10.8 


204 


18.3 


28.7, 109 


1+1 


7.0 


122 


18.3 


22.2, 84.0 




2+1 


n no 


13.0 


0.001 


36.9, 36.9 




1+0 


990 


2.62 10^ 


4.6710" 


49.3, 110 


DDC 


2+0 
1+1 


10.8 

12.5 


2.29 10^ 
1.89 10"* 


206 
103 


55.9, 122 
59.6, 212 




z+i 




ol / 


0.005 


126, 178 




1+0 


2.8 10^ 


7.45 10 


4.6710" 


83.0, 186 


DDC4^ 


2+0 


40.7 


8.64 10 


206 


109, 237 


1+1 


58.9 


8.59 10 


103 


127, 452 




2+1 


io.y 




0.005 


285, 403 




1+0 


2.8 10^ 


7.45 10 


1.0610" 


120, 269 


DDC4,r 


2+0 


40.7 


8.64 10 


10.6 


211, 791 


100 Hz 


1+1 


58.9 


8.59 10 


5.31 


242, 1116 




9-1-1 


13.9 


1.63 10^ 


0.00005 403, 403 




1+0 


13.010^ 


8.010'' 


1.0610" 


394, 880 


100 Hz 


2+0 


377 


1.8610^ 


10.6 


980, 3672 


Tfd = 5MeV 


1+1 


582 


1.9810^ 


5.31 


1160, 5354 


high Li^^ 


2+1 


171 


46.610^ 


0.00005 


2158, 2158 



TABLE I: Sensor effective volume for the Garching spectrum, 

^-av = 17.1 MeV, various detector configurations and detec- 
tion modes. The last entries, corresponding to DDCY^, are 
results for an enhanced DDC configuration with larger mod- 
ules and Q.E., for a T = 5 MeV Fermi-Dirac spectrum and a 
luminosity Ly^ normalized to 0.5 • 10^^ erg. 



o 




20 40 60 80 100 120 140 




-100 -50 50 100 150 200 

18 1| [ns] 

FIG. 6: Difference in hit arrival time. Distributions peaking 
from left to right correspond to the DDC4,r, DDC, DC and 
IC configurations. 



TITAND 49J project discussed in [T3] and designed for 
routine detection of supernova explosion at a rate of the 
order 1/yr. 



B. Energy spectrum 

The hit ratios (1 + 1)/(1 + 0) and (2 + 0)/(l + 0) are 
illustrated in Fig. [7j Both observables are sensitive to 
the average positron energy and do exhibit a 4.7%/MeV 
variation around the Garching energy for IC. All config- 
urations exhibit equally strong variations. In addition to 
their sensitivity to the spectral parameters of the emis- 
sion spectrum, observables based on ratios have the ad- 
vantage of canceling some of the involved systematics. 

For the benchmark Garching model. Fig. [8] provides 
the mean values for the neutrino spectrum, the positron 
spectrum and for IC, for the various coincidental modes. 
It is noted that the mean energy increases with the mode 
order. 

The resolution on the spectrum will depend on the 
accumulated statistics and can easily be calculated from 
the rig values presented in Table|l]and the distance. Con- 
sidering f^i {Nio 3> A^ii), the statistical uncertainty is 
approximately l/^/Ny^. Assuming a spectrum obeying 
Sfl model with a slight uncertainty on (E^), we can es- 




15 20 
<£e> [MeV] 

FIG. 7: Coincidental hit ratios (l + l)/(l + 0) (thin hues) and 
(2 + 0)/(l + 0) (think Hues) for DDC4^, DDC, DC and IC, 
and w.r.t. {E^+) (the same a was chosen for all {E^+)). 



timate A{E^) from the measurements, 



driydE{{E^))VWi 



(28) 



For example, the statistical error on the fraction 1 + 
1/1 + of the signal recorded in the IG configuration 
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conf. 


7 7 1 

kk 


T>Sf 


T->FD, T=5 McV 


fVFD, T=6.5 MeV 




1+0 


764 


1.0310'' 


1.32 10'' 


IC 


2+0 


7.4 


11.5 


16.4 


1+1 


1.67 


3.00 


4.83 




2+1 


0.12 


0.25 


0.45 




1+0 


990 


1.3310^ 


1.710^ 


DC 


2+0 


10.8 


16.8 


24.0 


1+1 


7.0 


12.2 


19.0 




2+1 


0.98 


1.90 


3.20 




1+0 


990 


1.3310^ 


1.710^ 


DDC 


2+0 


10.8 


16.8 


24.0 


1+1 


12.5 


20.6 


30.8 






AO 
Z.DZ 




; .z4 




1+0 


2.810^ 


3.810^ 


4.810^ 


DDC4,r 


2+0 


40.7 


63.0 


90.0 


1+1 


58.9 


93.2 


135 




2+1 


13.9 


23.4 


35.2 



TABLE II: Comparison of the sensor effective volume in m^ 
for Sf and Fermi-Dirac (T = 5 and T — 6.5 MeV) spectra for 
the various detector configurations and detection modes. 



for a benchmark SN assuming the Sf spectral shape and 
located at the galactic center is about 5% and does there- 
fore correspond to a resolution on the average energy of 
about an MeV (and even 0.2 MeV for DDC4^). 

This result is extremely promising and deserves a more 
sophisticated analysis beyond the scope of this paper in 
order to study the potential for characterizing the spec- 
tral shape of the neutrino flux, and possibly its full evolu- 
tion, i. e. the potential for obtaining the spectral neutrino 
light curve following a galactic SN. 



C. Neutrino flux directionality 

In the case of a galactic SN, it is desirable to locate 
the provenance of the SN neutrinos to inform the com- 
munity, which could immediately concentrate various in- 
struments toward this target to observe the delayed opti- 
cal burst. Moreover, it is likely that some supernovae are 
screened and never detected optically - pointing to the 
target in advance would result in a large multi-messenger 
sensitivity gain. 

The underlying idea relies on the combination of 

1) the fact that the average positron energy leading to co- 
incidental hits is higher than for single hits, as illustrated 
in Fig. [Sj which shows the detected positron spectra for 
the various modes. There is a gradual increase of the 
average energy of the modes with an increasing number 
of hits. 

2) an anisotropic scattering of the interaction with a 
slight preference for photons emitted forward if the higher 
energy tail of the neutrino emission can be selected. The 
average interaction {cos 9) is quite small around SN av- 
erage neutrino energies. The neutrino interaction is for- 
ward scattered above « 14 MeV, which represents be- 
tween 80 and 90% of the detected neutrinos depending 
on the coincidental detection mode). 



(cos 6') = 1% in the 1 + detection mode and rises to 
approx. 2% in the 1 + 1 detection mode with IC [38] (see 
the average energies in Fig.|8|. In order to exploit this 1% 
difference, the statistical uncertainty on fj^' must be kept 
smaller and requires therefore large statistics only avail- 
able below about 2 kpc with IC and 1 kpc with DC (and 
the azimuth of the incoming direction is not resolved). 
The situation improves sig nificantly with DDCY^ which 
will provide directional information for a SN exploding 
anywhere in the MW, moreover in a design with small 
inter-string spacing, the azimuthal direction will be con- 
strained as well (our ID framework must be extended in 
order to include inter-string coincidences). 
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Neutrino (thin, dotted) and 
produced (thin, dot-dashed), pro- 
double hit (dashed), 1 + 1 (dot- 
dashed) and 2 + 1 (dotted) signatures for IC. The resp. mean 
energies (in MeV) are: 
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positron spectra (Eq. |18[ ) 
ducing single hit (solid) 
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It is also clear that the time difference distribution 
of detected coincidental hits, an observable independent 
from f^x^ should also depend on the angle of provenance 
of the neutrino flux t?. 

A careful assessment and quantification of the resolu- 
tion by studying the impact of the incoming direction on 
the induced distortion on the hit time difference distri- 
bution and the various coincidental hit ratios requires 
Monte Carlo studies. Such a study should also take 
into account the elastic scattering channel. Although 
we pointed out above that the low cross section of this 
channel will limit its contribution to the overall detec- 
tion potential, the strong angular correlation and forward 
peaked cross section might result in a slight enhancement 
for the measurement of the incidence direction. 



V. CONCLUSIONS 

In order to confirm the simple analytical calculations 
for single hit effective detection volume and to precisely 
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study the new coincident hit detection methodology and 
the new observables to which it grants access, precise 
simulation should be performed. In particular, a precise 
assessment of the angular resolution in locating the SN 
from the asymmetric Cherenkov emission can be reliably 
studied as a function of distance and which is only qual- 
itatively stated in this paper. 

With such a tool, a dedicated work not focused only 
on IBD but implementing other neutrino interactions, 
e.g. with oxygen, would permit the precise assessment of 
the potential of the novel analysis method for oscillations 
and matter effect, hierarchy and deleptonization neutrino 
emission phase studies. 

A significant increase in potential would be reached 
with the imagined DDCj^ IceCube extension, with sen- 
sors combining several improvements such as a larger 
and uniform collection surface and an increased photo- 
detection efficiency (which is of specific interest as the 
distance reach scales linearly at first order with the Q.E. 
for 1-1-1 and 2-1-0 detection modes) . SN detection up to 
2 Mpc would be at hand with the foreseen potential boost 
of this detector. Such an extension, considerably cost ef- 
fective compared to similar projects, would be motivated 
by a broader experimental interest, vastly enhancing the 
detection potential of low energy neutrinos, > 1 GeV, 
connected to physics of neutrino oscillations, dark matter 
neutrino signatures, astrophysical neutrinos in particular 
from the southern hemisphere, etc [5T]. 

To reach about 2 Mton effective detection volume, thus 
granting access to SN at distance scales larger of about 
4 Mpc and guarantee routine SN detection, the DDCj^ 
configuration is likely too modest and the number of 
modules should rather be about lO'OOO. Also R&D and 
extensive simulations will be necessary in order to deter- 
mine the optimal module arrangement, the sensor tech- 
nology (e.g. the use of solid state detectors instead of 
PMT, see [ST] for a review), and other hardware related 
parameters. 

It is worth mentioning an interesting outcome of this 
analysis: the method described here can be readily ap- 
plied to existing and future IceCube data in the regular 
muon data acquisition mode: events of > 20 /zs duration 
are recorded at a rate of about > 2 kHz, i.e. representing 
a lifetime fraction of about 5%, without the necessity of a 
dedicated data stream. Numbers provided in this paper 
can readily be rescaled to obtain the reach in the base- 
line IceCube configuration and its DeepCore provided a 
corresponding rescaling of the distance reach (by a factor 
of about 4). Concerning a dedicated data stream, we do 
not discuss here the possible hardware implications for 
the realistic implementations of the methodology, but it 
may well already be implementable within the current 
IceCube data acquisition hardware. 

The methodology is natural for operating |52j and pro- 
jected [S3] water detectors, given a comparatively large 
noise rate (based on [5S], we found in [55] a discussion 
based on a similar approach). 

We finally note that accounting for spin-flavor conver- 



sion and MSW effect for neutrinos propagating outward 
the SN core [54] , the D^. flux is shifted to slightly higher 
average energy, improving the prospects of the coinciden- 
tal hit method. 
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